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Within generalized linear response theory, an expression for the dielectric function is derived 
which is consistent with standard approaches to the electrical dc conductivity. Explicit results are 
given for the first moment Born approximation. Some exact relations as well as the limiting be- 
haviour at small values of the wave number and frequency are investigated. 



I. INTRODUCTION 

The dielectric function e(k,u>) describing the response of a charged particle system to an external, time and space 
dependent electrical field is related to various phenomena such as electrical conductivity and optical absorption of 
light. In particular, it is an important quantity for plasma diagnostics, see, e.g., recent applications to determine the 
parameters of picosecond laser produced high-density plasmas [Q. However, the application of widely used simplified 
expressions for the dielectric function is questionable in the case of nonideal plasmas. 

As well known, the electrical dc conductivity of a charged particle system should be obtained as a limiting case of the 
dielectric function. However, at present both quantities are treated by different theories. A standard approach to the 
dc electrical conductivity is given by the Chapman-Enskog approach |2j. In dense plasmas, where many-particle effects 
are of importance, linear response theory has been worked out to relate the conductivity to equilibrium correlation 
functions which can be evaluated using the method of thermodynamic Green functions, see H. This way it is possible 
to derive results for the conductivity of partially ionized plasmas not only on the level of ordinary kinetic theory, but 
to include two-particle nonequilibrium correlations as well 

On the other hand, the dielectric function can also be expressed in terms of equilibrium correlation functions, but 
the systematic pcrturbativc treatment to include collision effects is difficult near the point k = 0, u> — 0, because 
an essential singularity arises in zeroth order. Different possibilities are known to go beyond the well-known RPA 
result. In the static limit, local field corrections have been discussed extensively p], and the dynamical behavior of 
the corrections to the RPA in the long-wavelength limit was investigated in the time-dependent mean field theory 
neglecting damping effects 0], see also for the strong coupling case. At arbitrary k, ui, approximations are made 
on the basis of sum rules for the lowest moments However, these approximations cannot give an unambiguous 
expression for e(k,to) in the entire k, lu space. 

We will give here a unified approach to the dielectric function as well as the dc conductivity, which is consistent 
with the Chapman-Enskog approach to the dc conductivity and which allows for a perturbation expansion also in 
the region of small k, lu. In the following Section II the method of generalized linear response ^ is presented which 
allows to find very general relations between a dissipative quantity and correlation functions describing the dynamical 
behaviour of fluctuations in equilibrium. A special expression for the dielectric function is given in Section III which 
is related to the use of the force-force correlation function in evaluating the conductivity. 

Different methods can be applied to evaluate equilibrium correlation functions for nonideal plasmas. We will use 
perturbation theory to evaluate thermodynamic Green functions Jl0| |. Results in Born approximation are given in 
Section IV. Using diagram techniques, partial summations can be performed as shown in Ref. ||. An alternative to 
evaluate equilibrium correlation function in strongly coupled plasmas is given by molecular dynamical simulations. It 
is expected that reliable results for the dielectric function for dense systems by quantum molecular dynamics will be 
available in the near future. Works in this direction are in progress but will not be discussed in this paper. 

To illustrate the general approach, explicit results for the dielectric function in lowest moment Born approximation 
are given for a Hydrogen plasma in Section V. A sum rule as well as the conductivity are discussed. The simple 
approximation considered here will be improved in a subsequent paper JTlJl , where a four- moment approach to the 
two-component plasma is investigated. 
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II. DIELECTRIC FUNCTION WITHIN GENERALIZED LINEAR RESPONSE THEORY 



We consider a charge-neutral plasma consisting of two components with masses m c and charges e c , where the index c 
denotes species (electron e, ion i) and spin, under the influence of an external potential U ext {r, t) — e % ^ kr ~ wt ^U ex t(k, w) 
+ c.c. The total Hamiltonian H tot (t) = H + H cxt (t) contains the system Hamiltonian 

H = E E v 4 c v + \ E V <*M c U 4+, d V c p (!) 

c.p cd.pp'q 

and the interaction with the external potential 

tfext(i) - U ext (k,w)e- iut J2 e - n P,-k + cc, (2) 

c.p 

where = h 2 p 2 /2m c denotes the kinetic energy, V c d(q) — e c ed/(e £lo? 2 ) the Coulomb interaction and Slo the 
normalization volume. Furthermore we introduced the Wigner transform of the single-particle density 

n p,k = ( n p,-k) = C p-k/2 C p+k/2- ( 3 ) 

Under the influence of the external potential, a time-dependent charge density 

7^ E ec <K,fc'>' + c - c - = 7T E e - W; w ) ei( ^" wt) + c - c - ( 4 ) 

"0 ^T, "0 — 

c,p,fc' c.p 

will be induced. Here, <5n^ fe , = n£ fe , - Tr |rip fe , poj denotes the deviation from equilibrium given by 

Pa = exp(-/3# + ^ cN ^ / Tr e M-PH + 1* 1 **) ■ ( 5 ) 

c c 

The average (...)* = Tr {.../>(£)} has to be performed with the nonequilibrium statistical operator p(t), which is 
derived in linear response with respect to the external potential in Appendix A. For homogeneous and isotropic 
systems, we find simple algebraic relations between the different modes (k,u>) of the external potential U ex t(k,u>) and 
the induced single-particle distribution 

5f c (p;k,co) = e iut (Sn^)* (6) 

which allow to introduce the dielectric function e(k,u>), the electrical conductivity a(k,u>), and the polarization 
function U(k,u>). From standard electrodynamics we have 

i 1 

e(k,u>) = H a(k,u>) = 1 — U(k,u>) , 

e uj e k z 

m^) = ^e c S m k,.) u ^ (7) 

with U e s(k,cj) = U ex t(k,cj)/e(k,ui). Using the equation of continuity 

k -> 

UJ ^2 s fc{p;k,u) = — E^ z S fc(P^ k ^) > (8) 

7Tl c 

p p 
where the z direction is parallel to k, k = ke Zl we can also express 

II(fc, w) = - ^- V — hp z $fc(p; k, u) 

uj U e s(k,w) 
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with the current density operator 

^E^-M- (10) 

c,p 

The main problem in evaluating the mean value of the current density (Jfc)*, Eq. (|l0|), is the determination of p(t). In 
linear response theory where the external potential is considered to be weak, the statistical operator p(t) can be found 
as shown in Appendix A. An important ingredient is that a set of relevant observables {B n } can be introduced whose 
mean values (B n ) characterize the nonequilibrium state of the system. The non-equilibrium statistical operator is 
constructed using a corresponding set of thermodynamic parameters (/> n (t). For weak perturbations, in linear response 
theory it is assumed that the <p n (t) are linear with respect to the external potential, and a set of generalized response 
equations is derived which allow to evaluate the response parameters 4> n (t) . The coefficients of these response equations 
are given in terms of equilibrium correlation functions which can be evaluated using the methods of quantum statistics. 

Solving this set of linear response equations by using Cramers rule, the response parameters can be eliminated. If 
the current density operator Jk can be represented as a superposition of the relevant observables B n , we find 

k 2 

u{k,uj) = i— pn 

LO 

with 

M 0n (k,Lj) = {J k \B n ) , M m0 (k,uj) = (B m ;J k ) , 

M mn {k, lj) = (B m ; [B n - iujB n ]) + (B m ; [B n ~ iuiB n ]) u+vri - J ^ u+lv (B m ; [B n - iu>Bn\)u+i V ■ (12) 
The equilibrium correlation functions are defined as 

i r p 

{A- B) = (B+; A+) = ~ dr Tr [A(-ihr)B+ p ] , 
P Jo 

poo 

(A;B) Z = / dte* zt (A(t);B) , (13) 
Jo 

with A(t) = exp(iHt/H) A exp(—iHt/h) and A = i[H,A] , furthermore we used the abbreviation 

j k = e- 1 (k,Lj) J k . (14) 

The correlation functions can be evaluated by standard many particle methods such as perturbation theory for 
thermodynamic Green functions. In this context the correlation functions containing Jk are obtained from irreducible 
diagrams to Green functions containing Jk, which do not disintegrate cutting only one interaction line. 

The expression ( [Tl|) for the polarization function is very general. Depending on the set of observables {B n }, different 
special cases are possible such as the Kubo formula or the Boltzmann equation to be discussed in the following section. 
It is also possible to include two-particle nonequilibrium correlations if an appropriate set of B n is chosen. We 
will work out here an approach to the dielectric function which is closely related to the Chapman-Enskog approach 
to the electrical conductivity. 



M 0n (k,Uj) 

M m0 {k,u) M mn (k,uj) 



\M mn (k,u)\ 



(11) 



III. MOMENT EXPANSION OF THE POLARIZATION FUNCTION 



Up to now, B n was not specified. It is an advantage of the approach given here that different levels of approximations 
can be constructed, depending on the use of different sets of B n . If no finite order perturbation expansion of the 
correlation functions is performed in evaluating the polarization function fliT| ) , ah these different approaches are exact 
and should give identical results. However, evaluating the correlation functions within perturbation theory, different 
results for the polarization function are expected using different sets of B n . As has been shown for the electrical 
conductivity results from finite order perturbation theory are the better the more relevant observables are 

considered. 

A simple example for a relevant observable B n characterizing the nonequilibrium state of the system is the current 
density (|l0|), 
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B n = J k . (15) 

During this paper, we will treat this approach in detail. The current density is related to the lowest moment of the 
distribution function. Possible extensions to more general sets of relevant observables are discussed at the end of this 
section. 



In the approach given by Eq. (15), we have 



with 



n M = — , (16) 



Mjj = —iuj (J k ; Jk) + (jk] Jk)u+vq — j[ (j fc . j k ) u+iv . (17) 

\Jk ) J k) u+ir] 



For the derivation we used the property 



(A;B) = ^Tr{[A,B+}po} (18) 



(for proving perform the integral in the definition ( |13 | )) so that (J k ; Jk) = jh§ Tr{[Jfc, J- k ]po} = 0. 
Applying integration by part ( |57| ) , the expression (|l6|) can be rewritten as 

jT/i ^_ ik 2 /3fl (JfcjJfc) (Jk\Jk)uj+iri n „N 

W { J kt J k) ~ V Wfe) J k)u+in 

Performing the limit 77 — )• , for finite values of the correlation function (Jk; Jk)^+ir) we obtain the simple result 

n(fc >w ) = ^ {Jk',Jk)*+iT, (20) 

LO 

which is also denoted as the Kubo formula for the polarization function. Similarly, the Kubo formula can also be 
obtained from more general sets of observables {£?„}. A direct derivation of the Kubo formula is obtained from 
Appendix A, Eq. (p2|), if the set of relevant observables B n is empty. Different approaches based on different sets of 
relevant observables B n are formally equivalent as long as no approximations in evaluating the correlation functions 
are performed. 

However, expressions ([Hj] ) and ( |20| ) are differently suited to perform perturbation expansions. For this we consider 
the static conductivity a — cr(0, 0) which follows from 



a(k, u) =i ^ n(fc,w) (21) 



in the limit k — > 0, uj — > . 
Comparing the Kubo formula 



with the result according to (fL6|), 



<T = {3n Q (J ; Jq) iv (22) 

(joj ^0) {Jo, Jq) £23) 

(•^0; Jo/i-n — (Jo; Jo)in (Jo; Jo) in (Jo; Jo)^ 

it is evident that perturbation theory cannot be applied to ( ^2|) because in zeroth order this expression is already 
diverging. In contrast, ( p3| ) allows for a perturbative expansion. For instance, in Born approximation the Faber - 
Ziman result for the electrical conductivity is obtained. The expression a^ 1 ~ (Jo; Jo)irj is also known as the force- 
force correlation function expression for the resistivity. More precisely, the resistivity should be given in terms of 
stochastic forces which are related to the second term in the denominator of Eq. (|23|), see also Eq. (|5|) in App. A. 
The applicability of correlation functions for the inverse transport coefficients has been widely discussed || . 

The approach to the dielectric function given in the present paper is based on the choice (|15| ) for the set of relevant 
observables and may be considered as the generalization of the force-force correlation function method for the electrical 
resistivity to the dielectric function. Possible extensions of the set of relevant observables have been investigated in 
evalua ting the dc conductivity in Ref. || and will be considered in evaluating the dielectric function in a forthcoming 
paper 
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IV. EVALUATION OF CORRELATION FUNCTIONS 



Within the generalized linear response approach, the polarization function is given in terms of correlation functions 
which, in general, are elements of matrices. Within a quantum statistical approach, the correlation functions are 
related to Green functions which can be evaluated by diagram techniques. This has been discussed in detail in the 
case of the static electrical conductivity H and will not be detailed here. Instead, we will consider only the lowest 
orders of perturbation theory (Born approximation). 

In the case considered here, the relevant observable Jk ( |To| ) is given by a single particle observable. The correlation 
functions occuring in ( [l6|) will contain the operators n c p k = Cp_ k / 2 c p+k/2 an< ^ k = ~ (^Pzk / m c ) n c p k + v° k , with 



u p,k 



d,p',q 



C p-k/2-q dp'+q d p ' C p+ k/2 C p-k/2 ^p'+q ^p' C p+k/2+q 



(24) 



To evaluate the correlation functions, we perform a perturbation expansion with respect to the interaction V, see 
App. B. In addition to the zeroth order terms, which reproduce the RPA result, we consider the Born approximation. 
Up to second order with respect to V we have 



m ^{J k ;J k \ 



i i v 

1 ^cd.pp 



pp' OS m c m d 



PzP 



ri-iuj+i-^p'^k rj—iui+i-^p^k (Jfc;Jfc) (0) 



(25) 



The evaluation of the correlation functions for the non-degenerate case is shown in App. B. We obtain the following 
expression 



n(fc,w) = -- 



l3J2 c e 2 c n c [l + z c D(z c )] 



1 - % 



2(2tt) 1 / 2 



Jfdper* 111^ + ^ )W(p) 



(26) 



with 



TO = ~ M — ~ — 



m; 



J2 C eg n c [l + z c D(z c ) 
Y, c e 2 c n c /rn c 



M. 



1/2 



1/2 



2L 



dc c 



e e D(z e 



-cp) + eiD{z el + 



-cp) 



(27) 



Here, z e 



2k B T> * c 



A(p) = (h 2 K 2 )/(4fi et k B Tp 2 ) + 1 , and 



D(z) 



dx 



x — z — irj 



— = i^e~ z [1 + Erf (w) 



(28) 



denotes the Dawson integral. Note that a statically screened potential was used in (53) to obtain a convergent collision 
integral, the screening parameter is given by k 2 = ^2 c e 2 n c / (eoksT) . From ( |2"6| ) it is immediately seen that the RPA 
result is obtained in the limit of vanishing interactions, W(p) = 0. 



V. RESULTS FOR HYDROGEN PLASMAS 



The expression ( |2q ) for the polarization function is simplified for a system consisting of protons and electrons, 
where = — e e , rij = n e , and mi/m e — 1836: 



e(k,w) = 1 + 



e k B Tk 2 



[2 + z e D(z e ) + Zl D(zi)} 



1 - i 



e4 n ^ 2(27r) y 2 r dp e V ( U *Z± + _A_ 

k*(4ne Q )* n (k B T)s/^ 2V) J dp \ A + 1 A + 1 



2 f M ■\ 1 ^ 2 f 1 ( Pm7 fm~ 

- p [2 + z e D(z e ) + ZiD(zi)] - — / dcc(D(z ei -J—cp)-D(z ei + J — 



cp) 



(29) 
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We first discuss the limiting case of small k. For k <C u)^Jm e / (2ksT) we use the expansion 



D(z) = iV^e~ z2 - - - _ ± . . . (30) 

z 2z 6 



so that after expanding also with respect to cp/ z e i we have 



e(0,a;) = l- pl - (31) 

0J Z + tUJ/T 



with = e 2 n/(eo/x e i) and 

_ _ (4ne )* {k B Tfl^ 3 



/ dp p e~ p (In 

vo 



A — 1 



A + 1 A+ V 



(32) 



e 4 n 4(2tt) 1 /2 

According to (j^), the dc conductivity 

CT(0,c;-fO) = w^ 1 eoT (33) 

is obtained, what coincides with the Faber-Ziman formula at finite temperatures jjj. 

On the other hand, in the limiting case of small u> we use for oj <C y // 2kBT/rriik the expansion 

£>(z) = iv^e -22 - 2z + \z 5 ± . . . (34) 

and obtain 



with 



lim lim e(k,u) = 1 + - = ( l + i — J^- ) (35) 



(47reo) an (fc S T)5/2 / p A + 1 A + 1 W 



Here, in evaluating the last expression of (|29|), also z e i + m e / rriicp is considered as a small quantity, whereas 
Zpa — \/ rrii I m e cp is large in the region of relevant p. For small values k < ^/2fcaT/(7rroi)2/<j, the second term in the 



numerator of ( p5| ) can be neglected, and the diffusion type form of e(k, to) is obtained, see |12 |. 

As an example, a dense plasma is considered with parameter values T = 50 eV and n e — 3.2 10 23 cm' 3 . Such 
parameter values have been reported recently in laser produced high-density plasmas by Sauerbrey et al. , see 0] . We 
will use Rydberg units so that T — 3.68 in Ryd and n e = 0.0474 in a^ 3 . At these parameter values, the plasma 
frequency is obtained as cj p j = 1.54, and the screening parameter as k = 0.805. 

First we discuss the dependence of the dielectric function on frequency for different values of k, see Figs. 1-4. For 
large values of k our result for the dielectric function coincides with the RPA result. At decreasing k strong deviations 
are observed. Both the RPA expression as well as the expression (^9|) for the dielectric function fulfill important 
relations such as the Kramers-Kronig relation and the condition of total screening. The validity of the sum rule 



u lme(k, u) dui = ■^pi (37) 



is checked by numerical integration. The RPA result coincides with the exact value ut^n/2 = 3.74 to be compared 
with expression (§|) which gives 3.74 at k = 1, 3.75 at k = 0.1, 3.71 at k = 0.01 and 3.74 at k = 0.001. The small 
deviations are possibly due to numerical accuracy. 

To investigate the behavior at small k, we give a log-log plot of lme(k, w) as function of w for different values k in 
Fig. 5. For uo > ^/2fcgT/m e fc = 3.84fc the Drude-like behaviour (|3l]) is clearly seen, with r = 8.36. 

Considering the limit of small uj, a log-log plot of Ime(fc, lu) as function of k for different values u> is shown in Fig. 6. 



The diffusion behavior (j3|) occurs for k < v / 2k B T/ (nm,) = 0.00732 at k > ^/m,/ (2k B T)ui = 11.17a; with d = 13.8 . 
Altogether the numerical evaluation of the general expression (^9[) for the dielectric function confirms the validity of 
the simple limiting formulae (|3l]) and (p5[). 

In this paper we have focussed the discussion only to the properties of e(k,u>). Related quantities such as e~ 1 (k,uj) 
will be investigated in a forthcoming paper p"l[ . The parameter values for density and temperature can be extended 
to other nondegenerate plasmas like ordinary laboratory plasmas or the solar plasma. This has been done with results 
showing the same qualitative behavior of the expression (^9|) in comparison with the RPA expression, but at shifted 
values of k and ui. 
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VI. CONCLUSIONS 



An expression for the dielectric function of Coulomb systems is derived which is consistent with the Chapman- 
Enskog approach to the dc conductivity. For a two-component plasma, explicit calculations have been performed in 
the lowest moment approach. In Born approximation, expressions are given which allow the determination of e(fe,oj) 
in an analytical way. It is shown that general relations such as sum rules are fulfilled as well as the dc conductivity 
is obtained in the form of the Ziman-Faber result. 

We performed exploratory calculations to illustrate how the generalized linear response approach works. Obviously 
an improvement of the results can be obtained if i) the Born approximation is improved including higher order of 
perturbation theory, ii) higher moments of the single-particle distribution are taken into account. Both points have 
been discussed for the limiting case of the dc conductivity || , where a virial expansion of the inverse conductivity 
was given. 



A four moment approach will be presented in a subsequent paper 1 1 where also the comparison with the Kubo 
approach and computer simulations are discussed. Within the approach given here it is also possible to treat the 
degenerate case. Work in this direction is in progress. 
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APPENDIX A: GENERALIZED LINEAR RESPONSE THEORY 

To construct the nonequilibrium statistical operator p(t) we use the density matrix approach |^,^3| . Characterizing 
the nonequilibrium state of the system by the mean values (A„) f of a set of relevant observables the generalized 

Gibbs state 

p Tel (t) = e~ s W kB , (38) 
-S(t) = $(t)+ Y, c*n(t)A n (39) 



where 



is the entropy operator and 



*(*) =lnTrexp|-^a„(i)A„| 



(40) 



is the Massieu-Planck function, follows from the maximum of the entropy 

(S(t)>* = -k B Ti{ Prcl (t) lnp re i(t)} (41) 

at given mean values 

Tr{A n p lel (t)} = (A n y . (42) 

The thermodynamic parameters (Lagrange multipliers) a n (t) are determined by the self-consistency conditions ( |42| ) 
and will be evaluated within linear response theory below. 

The relevant statistical operator ( |38| ) does not solve the von Neumann equation, but it can serve to formulate the 
correct boundary conditions to obtain the retarded solution of the von Neumann equation. Using Abel's theorem, 
the nonequilibrium statistical operator Q is found with the help of the time evolution operator U{t, t'), 

ih^- t U(t,t')=H tot (t)U(t,t'); U{t',t') = l, (43) 

as 



7 



p{t) =r, dt 1 e-**-*'' U(t, t') p Tel (t') U(t', t) , (44) 



where the limit r\ — * has to be taken after the thermodynamic limit. Partial integration of (|44| ) gives 

P(t) =Prel(i) + /W(*) (45) 

with 



/Wi(t) = - ^ dt' e-"(*-^l7(*,f) | ~ [H tot (t') , Pxel (t')} + JL Prel (t')\ tf(f >t ). 



(46) 



The self-consistency conditions (^2|) which determine the Lagrange multipliers take the form 

(*)} = . (47) 

For a weak external field U cx t, the system remains near thermal equilibrium described by po (|^), so that p(t) (EE) 
can be expanded up to the first order with respect to U CX f For this we specify the set of relevant observables {A n } as 
{H, N c , B n (r)} (note that summation over n in ( (39| ) also means integration over r) and the corresponding Lagrange 
parameters {a n } as {ft, — /3/i c , —@<l>n('i f ,t)}, 

^-S(t) = m)+j3H-(3j2»cN c -!3j2 J d 3 r n (r,t) B n (r) . (48) 

c n 

Expanding the nonequilibrium statistical operator up to first order with respect to f7 ext and <fi n (r , i) it is convenient 
to use the Fourier representationR so that 



f d 3 r 0„(r, t) B n {r) = «) e"*" S+ + c.c. (49) 
with 

0„(r, i) = e*^— *) B (£, w) , B„ = J d\ B n (r) e"'^ . (50) 
The contributions to (EJ|) are 

PreiW -Po + e" Mt / frVB+(ffir)^(M)A) + cx. (51) 
Jo 



and, applying the Kubo identity 

[A,po] = / dre-™ [H,A] e™ p 
Jo 



(52) 



we find 

rt 



- ^ (*' - t + ihr) - iuB+ (t' -t + ihr) (f> n (k, w) \ 

n ) 



Po + c.c. (53) 



l In general we have (f> n (r,t) = Ylh> I ^ e' (fcV ^' l) 4> n (k' ,lu') and B n> y = J d s rB n (r)e lk? . The selfconsistency equations 
([I?]) must be fulfilled for any time t so that u' — u> follows. Furthermore, the equilibrium correlation functions Tr(j4fc_B^",po) do 
not vanish only if k! = k so that [Ak {r) — iu)) \ B k i ) ~ &kk> ■ The well-known property of linear response that only such fluctuations 
are induced where the wave vector and frequency coincide with the external potential is a consequence of homogeneity in space 
and time. 
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Inserting this result in the self-consistency conditions ( pLT] ) we get the response equations 

— (B m \ A) U j+i n U c s(k,uj) = (B m ;C) u+ i v 
with the correlation functions defined by (p"3|), 

A = ^e c n pk = ikn J k , 



and 



C = 



iu>B n 



(54) 
(55) 

(56) 



To make the relation between the response equations ( p4| ) and the Boltzmann equation more closely, see we 
introduce the 'stochastic' part of forces applying partial integrations 



{A;B) 2 

so that (|54|) can be rewritten as 

-iktto \B m \ J k ^j U cS (k, oS) 

= {B m \ C) + ( Br, 



(A;B) + (A;B) 2 



(A-B)-(A;B) 



(57) 



\^rn\ Jk) ~t~ {B m , Jk) uj-\-ir} {-E$mi Jk) cj+ir; 



(-^m; Jk) uj-\-iri 



B,, 



{Bfri! Jk) u)A-i'q 

q _ (B 7n '-> C)cj+ir] j 
{Bjn\ Jk) uj-\-ir} 



{B mi C)uj-\-ir/ 



{B m ^ Jk) LJ-\-ir] 

Then, we find the following form for the response equations 

- ikflo M m0 t/eff(fc,a;) = ^ M rnn (f) n (k,u) 



with 



and 



M, 



in — (yB m \ j k 



M mn = (B m ; [B n - iujB n ] 



r> _ {Bin] Jk)u+in j-, 

™ / D . T.\ n 
V^m; J k j uj+in 



; [B n - iujB n 



uj-\-ir) 



(58) 



(59) 



(60) 



(61) 



The system of equations ([39|) can be solved applying Cramers rule. Then, the response parameters are represented 
as a ratio of two determinants. 

With the solutions 4> n the explicit form of p(t) is known, and we can evaluate mean values of arbitrary observables. 
In^particular, we are interested in the evaluation of (Jk) 1 exp(iwi) to calculate the polarization function (|^) using ([5]]), 

(J k y e iut = f3j2 Bn) - (J^ [B n - iuB n ]) u+ir ,} 4> n {k,uj) 



ikttoP (Jk;Jk)u+in U c ft(k,u)) 



(62) 



If Jk can be represented by a linear combination of the relevant observables {B n }, we can directly use the selfconsis- 
tency conditions (M3) and have 



(J k y e iut =Tr[J fc Prcl (t)} e 1 ' 



(63) 



Comparing with ( p2|) we see that the remaining terms on the rhs of (|62|) compensate due to the response equations 
(|59|). After expanding p ro \(t) up to first order in (j> n (k,uj), Eq. (pl|), we have 



(Jk) 1 e*" 1 =/3^(J fc ;B„) 4> n {k,i 



(64) 



Inserting the solutions for 4> n in the form of determinants, we get the same result as obtained if we expand the 
numerator determinant (11) with respect to its first row. 
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APPENDIX B: EVALUATION OF THE COLLISION TERM IN BORN APPROXIMATION 



Let us first consider the lowest order of perturbation theory where we have for the correlation functions 

( n p,ki n p' ,k) — fp,k <W $cd 

( n i,k^ n p',k)^v = (r)-iv + i ^Pz k / m c)' 1 f^ k S pp' 6 cd , (65) 

where 

fp,k = m 2 Pz k/m c y\f^ k/2 - f£ +k/2 ) (66) 

Notice that lim^o/pfe = fp = {exp[/3(i?p — /i c )] + 1} _1 . In the classical limit where the Fermi function can be 
replaced by the Maxwell distribution, we have in lowest order with respect to the Coulomb interaction 

(J fc ;J fc ) (0) = %^E— ( 67 ) 



with Zc = l and 



Furthermore we have 



{J* J^X, = p ^ E e ' n ° t 1 + *°D{z c )] (68) 



D(z) = ± fe- 2 dx . . (69) 
V 7 *" J-oo x-z-ir) 



(A; *>&, = " ^ E I «• " £ £ E «! - 1' + ^Wl 

c c 
= — (Jki Jk)u+iTi j (70) 



(A; A-)! !,,, = -iw ^ ^- n c - i — e 2 c n c [1 + z c D{z c )\ , (71) 

so that from Eq. (jl^) the random phase approximation (RPA) 

(*, u) = -0^ e\ n c [1 + z c D(zc)} (72) 

c 

is obtained. 

After we have considered the collisionless plasma, we will now treat the general case of an interacting system 
where the correlation functions have to be evaluated with the full Hamiltonian ([!]). The evaluation of equilibrium 
correlation functions for an interacting many-fermion system can be performed within perturbation theory such as a 
Green function approach, and many-particle effects can be treated in a systematic way. We will give here the lowest 
order contribution with respect to the screened Coulomb interaction (Born approximation), a systematic treatment 
of higher orders can be done as indicated in || for the case of static conductivity. 

In the numerator of (|T6|), the higher order expansion for (J&; J&) lead to the replacement of the occupation numbers 
fp for the free fermion gas by the occupation numbers in an interacting fermion gas. This corrections in Born 
approximation can be given as shift of the single-particle energies and can be replaced by a shift of the chemical 
potential. 

We will investigate here the collision terms where the Born approximation leads to essential contributions. For this 
we use the relations (proof by partial integration (p7|)) 

KU; u p',fc)^+ t r, = (rt-iaj + ihpzk/mc)- 1 [{n c pk ;v^, k ) + (i£ fe ; V*, )k ) u+iv ] , (73) 
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{Vp,kl n t',k)u+in = (V + ihp' z k/m d ) 1 [(*£ fc ; r$ ifc ) - (t£ fc ; Up/ )fc ) w +t,] 
and find considering only the interaction in the collision terms 



(74) 



M 



33- 



(Jkj Jk) 2 
{J k 1 J k) uj+irj 

-1 



£ 



e c ed 



cd.pp' 
(Jk', Jk) 



f2p m c m d 



f j c d \ 
PzP z \ v p ,k> V p',k/u-\-i7) 



{J k j Jk)u-\-i7] 



l 



rj — ilu + ihp' z k/md rj — ito + iTip z k/m c 



(75) 



In the Born approximation for the frequency and wave vector dependent collision term we take the evolution 
operator due to the noninteracting part H° of the Hamiltonian (|l|) so that the correlation functions are immediately 
evaluated using Wick's theorem. Dropping single-particle exchange terms what can be justified for the Coulomb 
interaction in the low-density limit, we find 



(«p,fc(f7-^);up/, fe ) = £ 

c'p" q 

X {/p+fc/2-g(l - fp-k/2)3( E p+k/2~q 



-v ccl {q)f c pll+q {i-f;») 



TPC T?C 771C 

%"+q ~ ^p-k/2 ~ %" 



huj) 



x [V cc i{-q)Scd(5pi,p- q - Sp'.p) + V C 'c(-k + q)8c>d{5p',p"-k/2+q - V,p"-Hr/2)] 

~fp+k/2^ ~ fp-k/2+q)3( E p+k/2 + E p"+q ~ E p-k/2+q ~ E p" ~ ^ U ) 

x [V cc >(-q)6 c d{S P >, P - 8p>,p+q) + V c * c (-k + q)5 dd (8 p ,^„_ k / 2+q - 5 p >y, +q/2 )] } . (76) 

We evaluate the matrix element Mjj, Eq. ( f75|) in Born approximation to obtain the polarization function H(k,uj), 
Eq. (|l|). Using (|l7|), @, @, (|7§) we introduce 



i? = 



(Jk;Jk) { Z v 1 B " Y, c eln c [l + z c D{z c )) 



(77) 



and find the perturbation expansion Mjj = Mj°J + Mjj , where 



Mf} =R(J k ;J k ) 



(0) 



(78) 



m jj = 752 £ -7^ P*I* -1 + R 

"° cd,pl mcmd 



Tj — iuo -\- iHp z k/m, c r] — iu + ihl z k/md 
Evaluating the correlation functions (v'i k ; vf^u+i^ in Born approximation J76]), we have for small k, i 



(79) 



M j) - 2 § Q T, 5{E; +q + E\_ q -El- E\) qz (j± - ^) 



Ipq 

^ Pz + ^l z j - 2/7 

m e m, 



ihkp z /m e — iuj + i] m e ifikl z jm,i — iuo + rj m, 



(80) 



The further evaluation is done with introducing total and relative momenta P = p + 1, p 1 = (rriip— m e l) /M e i, p" = 
p 1 + q, M ei = m e + mi, fjQ 1 = m^ 1 + to" 1 so that 



Mj 1 ] = 2 



x e 



h-K e 2 e 2 ( 2irh 



2 \ 3/2 



rii 



Q ep \m e k B T 

h 2 p 2 k 2 p' 2 ,„ ,,„ 

2M ci k B T e 2 Mei fe fl T S(p — p ) 



2 \ 3/2 



1 2/^ e 



k B T) (2tt)9 h 2 



d 6 P I d 3 p' / d 3 p" 



1 



x P' z ~2R 



nip m, 



({f -p") 2 +K 2 ) 2 



{P'z ~ P Z ) 



TO e TO, 



ih 2 k 2 \ P 7 + M^v' - 4 



m, I z Uk 



(81) 



11 



Furthermore we introduce dimensionless variables HP(2M el k B P) 1/2 , hp 1 ' (2fjL ei k B T) 1/2 , A= (fi 2 k 2 ) / ' (^ ei k B Tp 12 )^ 
1 and spherical coordinates p' = {p'(l — c 2 ) 1 / 2 , 0,p'c}, p" = {p"(l — z 2 ) 1 ! 2 cos</>,p"(l — z 2 ) 1 / 2 sm(j),p"z} and perform 
the integral over <j) according to 



2tt- 



A — cz 



[\-cz- VT^VT^cos^] 2 (A 2 — 1 + c 2 — 2Acz + z 2 )3/ 2 



so that 



My J = —n e rii 



fl e 2 \{2^fk B P) 2 J p 



1 f°° 1 



7 dp'e- p / dc / dz 



r 3/2 



d 3 Pe~ p 



X — cz 



(A 2 - 1 + c 2 - 2Acz + z 2 ) 3 / 2 



(z-c) {p' 2 c — - — 



m e rrii 



iRp' — - — 



cj / M P . 



m e TOi / k B Tk 2 y /i e 



~ PC ~ k V 2fc B T ^ \ miP C k V 2fc B T 



Now, the integrals over z and P can be performed. Using 



(82) 



(83) 



. A — cz 

7 

! (A 2 - 1 + c 2 - 2Acz + z 2 ) 3 /: 



(z — c) = c I In 



we finally find 

r (i) 



1/2 



1 



2 / ( , ' 

3 "" \ m R m 



e 2 V2fc B r; 47r 3 / 2 7 

6 fi c 



dp e p I In 



A- 1 



A- 1 



A+l A+l 



'' 1 ' ,m e mi) k B Tk 2 v /j, ei J_ x 



<M F , 



A+l A+l 
i 



e e L>(z ei - » / — cp) + e l D{z ei + ./ — cp) 
m R V m, 



(84) 



(85) 



with z P i = ¥ 



M e 



el k V 2k B T 



Together with (|78|), (|67|), this result can be inserted in expression ( |l6| ) to evaluate II(fc, w). 
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Figure captions: 

Fig.l: e(k,uj) as function of u (in Ryd/?i) at k = 1/as for a hydrogen plasma, n e = 3.2 10 23 cm~ 3 , T = 50 eV. 
a: Re e, b: Im e. 

broken line: RPA, full line: first moment Born approximation. 
Fig. 2: The same as Fig.l for k = OA/as- 
Fig. 3: The same as Fig.l for k = 0.01/as. 
Fig.4: The same as Fig.l for k = 0.001/a B . 
Fig. 5: Im e(k,u>) as function of uj for different k. 
Fig.6: Im e(k,w) as function of k for oj = 0.000001 Ryd/h. 
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